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implantation and annealing. The NV distribution depends strongly on the native nitrogen concen- 
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tum information and magnetometer devices where NV centers must be created in close proximity 
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I. INTRODUCTION 

Nitrogen-vacancy (NV) centers located deep within 
a diamond lattice appear promising as solid-state spin 
qubits since they combine optical initialization and read- 
out capabilitie d * 2 * 3 *, long electron spin coherence life- 
times (approaching 1ms at room temperature^), an d the 
ability to control coupling to individual nuclear spins 5 6 7 . 
While these capabilities were initially demonstrated in 
high-purity natural diamond samples, similar results 
have rec ent ly been achieved in high-purity CVD-grown 
diamond^!. There is now much interest in fabricating 
spin-based devices in diamond, with potential applica- 
tio ns in quantum communic ation^!, qu antum computa- 
tior l 11 ' 12 1 , and magnetometr} 113 | 14 | 15 | 16 [ In all of these 
applications it is necessary or at least advantageous to 
couple NV centers to optical structures such as waveg- 
uides™! and microresonatorsSMniIimi23l2lll j to enable 

communication between distant qubits or to allow effi- 
cient extraction of emitted photons. Therefore, a reliable 
method is needed to create NV centers with good spectral 
properties in close proximity (< 100 nm) to a diamond 
surface. In addition, the charge state of the NV center 
must be controlled; all of the results mentioned above 
were demonstrated using NV centers in their negatively- 
charged state, NV~. 

It has been well established that NV centers can be 
created by ion implant ation or irradiation followed by 
annealing j 26 | 27 | 28 | 29 | 30 | 31 [ Ion implantation is used to cre- 
ate damage (carbon vacancies) in the diamond lattice, 
an d it can a lso be used to introduce nitrogen impuri- 
^ e j20 | 32 | 3 3i34|^ The choices of particle and acceleration 
voltage determine the initial vertical distribution of va- 
cancies. Annealing is performed at temperatures (typi- 
cally 600 — 1000°C) where vacancies become mobile and 
can combine with nitrogen impurities to form NV centers. 
Vacancies can also recombine with self-interstitials^ 7 ! 35 * 36 ! 
(implantation produces one interstitial carbon for each 
vacancy), can join together to form extended defects^, 
or can become trapped at dislocations or at the surface^. 
While numerous annealing experiments have been per- 



formed, and vacancy migration over long distances in the 
horizontal direction has been observed at 1400° C anneal- 
ing temperature^, to our knowledge none of these studies 
has yet clearly answered the question of how far vacan- 
cies diffuse vertically from the initially damaged layer at 
temperatures optimized to create NV centers. This is of 
vital importance for positioning NV centers close to an 
optical structure. 

Here, we present etching experiments that reveal the 
vertical distribution of optically active NV centers pro- 
duced using ion implantation and annealing to convert 
native nitrogen into NV centers. The results indicate that 
vacancies can diffuse several hundred nanometers to form 
NV centers at the annealing temperatures used. Differ- 
ent behavior was observed for high-nitrogen (w 100 ppm) 
and low-nitrogen (ss 1 ppm) diamond. In low-nitrogen di- 
amond, spectroscopy results show an unexpected depen- 
dence of the NV~ signal on etch depth, and we interpret 
this behavior as an electronic depletion effect. 



II. EXPERIMENT 

The first sample used in the experiments was a 3 x 3 x 
0.5 mm, (lOO)-oriented CVD diamond (Element 6) with 
a specified nitrogen concentration below 1 ppm, show- 
ing fairly uniform NV concentration before processing. 
The sample was implanted with 200 keV Ga ions (Core 
Systems) with a dose of 3 x 10 n cm~ 2 . A Monte-Carlo 
simulation using SRIM 40 predicts that this implantation 
energy produces vacancies to a depth of approximately 
100 nm, as shown in Fig. [I] During implantation, the 
CVD sample was covered with a TEM grid to allow 
measurement of the photoluminescence (PL) intensity 
contrast between implanted and non-implanted surfaces. 
This technique was necessary because the CVD sample 
has a high background NV concentration relative to its 
nitrogen content. After implantation, the sample was 
annealed at 925° C in an H^/Ar forming gas for 3 hours. 

The second sample was a 3 x 3 x 0.5 mm, (100) HPHT 
diamond (Sumitomo) with a specified nitrogen concen- 
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tration of 30 — 100 ppm. This sample has growth sectors 
with widely varying impurity concentrations. In most 
sectors, the width of the NV~ zcro-phonon line (ZPL) 
measured at low temperature is broad (~ 1 nm) , suggest- 
ing a high nitrogen concentration, and the NV lumines- 
cence achievable by implantation and annealing in these 
sectors is very high. However, as we have previously ob- 
served in Sumitomo HPHT samples^, there exist partic- 
ular sectors, with boundaries aligned to the sample edges, 
that have properties consistent with much lower nitrogen 
concentration, including a narrow ZPL lincwidth (as low 
as ~ 0.02 nm) at low temperature and an inability to 
produce strong NV luminescence with high implantation 
doses. The spectral properties of this low-nitrogen sector 
in the HPHT sample appear more similar to those of the 
CVD sample described above. The HPHT sample was 
implanted with 200 keV Ga ions (Core Systems) with a 
dose of 3 x 10 12 cm~ 2 . A TEM grid was not used on this 
sample, and was not required because the HPHT sample 
has a very low background NV concentration relative to 
its nitrogen content. After implantation, the sample was 
annealed at 925° C in an H 2 /Ar forming gas for 3 hours. 
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FIG. 1: Monte-Carlo simulations (SRIM) of the carbon va- 
cancy (black) and Ga (blue) profiles produced by 200 keV Ga + 
ion implantation. 
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FIG. 2: Fluorescence images of the samples obtained by scan- 
ning confocal microscopy before etching: (a) Element 6 CVD 
sample, showing implantation grid pattern; (b) Sumitomo 
HPHT sample, showing growth sector boundaries. Scan size: 
178 x 178 fjm 



The samples were etched using an Oxford inductively 
coupled plasma reactive ion etching (ICP-RIE) system 
with an O^/Ar recipe^. Before etching, protective masks 
(SiC>2 or SiN) were patterned on portions of the samples 
to allow measurement of the etch depth and for observa- 
tion of the PL contrast between etched and non-etched 
regions. After each trial, the etch depth was measured 
using a profilometer. Since the step height seen by the 
profilometer includes both the diamond etch depth and 
the mask thickness, for some trials we also measured the 
etch depth by electron microscopy. We measured the op- 
tical properties of the NV centers using a scanning confo- 
cal microscopy setup. For excitation, a 532 nm laser was 
focused through a 0.6 numerical-aperture microscope ob- 
jective to a ~ 0.5 [xia spot, and galvanometer mirrors or 
a piezo stage were used to scan the position of this spot 
on the sample. The resulting PL was collected through 
the same objective and imaged onto a pinhole for spatial 
filtering, and detected using a cooled-CCD spectrometer 
or an avalanche photodiode with bandpass filters select- 
ing wavelengths from approximately 647 — 804 nm, for 
detection of the NV~ phonon sidebands. The depth of 
focus was such that when the signal from a thin sheet 
of NV centers was measured as a function of the height 
of the microscope objective, a symmetric peak was ob- 
served with a full width at half maximum of 3.8 [Mm, cor- 
responding to 9 [Mm in diamond. The measurements were 
performed at room temperature. For both samples, im- 
ages of the PL intensity were obtained at the same lo- 
cations after each trial. For the HPHT sample, spectra 
were also obtained at multiple locations, and from these 
spectra the intensities of the NV° and NV - zero-phonon 
lines at 575 nm and 638 nm, respectively, were individu- 
ally determined. 

Fig.[2]shows PL images of the two samples before etch- 
ing. For the CVD sample, the implantation pattern is 
visible, combined with ribbon-like features that must be 
caused by a non-uniform nitrogen concentration present 
in the sample before processing. For the HPHT sample, 
the image shows a sharp boundary between growth sec- 
tors with different nitrogen concentrations. The region 
that appears dark is the low-nitrogen sector mentioned 
above. From images similar to these, the total photolu- 
minescence intensity in the 647 — 804 nm band was mea- 
sured as a function of etch depth for the two samples. 
The results, normalized by excitation power, are plot- 
ted in Fig. [3J For the CVD sample, the quantity plot- 
ted is the average intensity measured in the implanted 
squares minus the average intensity measured outside of 
the squares. In this way we can separate the PL pro- 
duced by NV centers created through ion implantation 
from the background NV PL. For the HPHT sample, the 
initial NV concentration was very small, and only the to- 
tal PL intensity is plotted in the figure, measured in both 
the high-nitrogen and low-nitrogen sectors. Since the op- 
tical detection efficiency in our setup is approximately 
constant over 1 [im of depth, the plotted quantities rep- 
resent the total PL intensity from all of the remaining 
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NV centers that were created through ion implantation 
and annealing. 

Before presenting a detailed analysis, we note several 
important features of the data in Fig. [3j First, the be- 
havior of the high-nitrogen sector of the HPHT sample 
for the first 300 nm of etching is quite different from 
that observed in the low-nitrogen sector of the HPHT 
sample and in the CVD sample. Considering just the 
HPHT sample, before etching, the total PL intensity 
is 44 times higher in the high-nitrogen sector than in 
the low-nitrogen sector, even through the implantation 
and annealing conditions were identical. As the sample 
is etched, the PL intensity in the high-nitrogen sector 
rapidly decreases, indicating that most of the NV cen- 
ters exist within the first 200 nm, while the PL intensity 
in the low-nitrogen sector is approximately constant, sug- 
gesting that the first 200 nm of material is devoid of opti- 
cally active NV centers. This is a strong indication that, 
for the high-nitrogen sector, the NV concentration pro- 
duced through implantation and annealing was limited 
mainly by available vacancies, while for the low-nitrogen 
sector the NV concentration was limited by the available 
nitrogen. However, at an etch depth of 400 nm, the PL 
intensities in the two sectors are almost equal (they differ 
by factor of 1.7), and at this point the NV concentration 
is probably limited by the available vacancies in both sec- 
tors. As the etching progresses further, the weak back- 
ground NV signals become important. Since the high- 
nitrogen sector has a higher background NV concentra- 
tion, the PL intensities in the two sectors again diverge 
for the largest etch depths. 



III. NV CONCENTRATION ESTIMATE 

The PL measurement integrates contributions from all 
of the remaining NV centers produced by ion implanta- 
tion and annealing, since the depth of focus of our opti- 
cal setup is much larger than 1 /im in diamond, as dis- 
cussed above. To a first approximation, the NV concen- 
tration is therefore proportional to the first derivative of 
the PL intensity with respect to etch depth (. However, 
there are two effects that complicate the analysis: spon- 
taneous emission modification by the nearby diamond-air 
interface, and electronic depletion. Spontaneous emission 
modification by a dielectric interface is well understood 
and can be included in the data analysis if we make an 
assumption about the radiative quantum efficiency of the 
NV center. In this section, we describe how to include 
spontaneous emission modification and obtain density es- 
timates that are corrected for this effect. Electronic de- 
pletion is more complicated, but appears to be important 
only for the low-nitrogen material. Electronic depletion 
is discussed and modeled in the following sections. 

The spontaneous emission rate of a quantum dipole is 
proportional to the total power radiated by an equivalent 
classical dipole. In the presence of a dielectric interface, 
the radiated power includes an interference term due to 



radiation reflecting off of the interface and performing 
work on the dipole. We calculated the ratio between the 
total radiated power, including the dielectric interface, 
and that in bulk diamond, using a formula derived in 
Ref. 1431 This ratio is equal to the spontaneous emission 
rate enhancement factor F p (z) (z is the depth below the 
surface) plotted in Fig [4j For the NV center, we must 
include two dipoles that are orthogonal to the NV axis 
and to each other. The NV axis can be oriented along 
any of the four (111) crystal directions, and in our sam- 
ples the surface is (100). For this geometry, the radiated 
power averaged over the two dipoles (regardless of how 
they are chosen) gives the same result as an isotropic 
mixture of dipole moments. For the result in Fig. |3Ja), 
we have also averaged over 650 — 750 nm emission wave- 
lengths, to take into account the broad NV - spectrum 
that is dominated by phonon sidebands. The sponta- 
neous emission enhancement (or suppression if F p < 1) 
is due entirely to the reflected wave, and only the ra- 
diated power into the substrate is modified. The radi- 
ated power into air is unchanged, but as F p increases, 
the spontaneous emission lifetime of the NV center de- 
creases, and thus the probability per excitation event for 
a photon to be emitted into air decreases. In the limit of 
weak excitation at 532 nm, the rate of excitation events 
does not depend on this spontaneous emission modifi- 
cation effect, and depends only on the driving electric 
field intensity, the dipole moment of the optical transi- 
tion, and the spectral width of the phonon sidebands. 
Thus, in the weak excitation limit, assuming a radiative 
efficiency of unity, the collection efficiency is inversely 
proportional to F p . We have defined a theoretical, rela- 
tive collection efficiency 7] — l/F p , plotted in Fig. [4|a). 
We have ignored one subtlety by using a single r\ based 
on the average radiated power from the two mutually or- 
thogonal dipoles of the NV center. If there were no relax- 
ation between excited states, then the two dipoles would 
each have a different spontaneous emission lifetime, and 
we would have to use r\ = 

( F p,i + F p,l)/ 2 - However, at 
room temperature, where rapid thermalization between 
the two orbital excited states may be expected, we be- 
lieve it is more appropriate to average F p before taking 
the inverse. 

Taking spontaneous emission modification into ac- 
count, we expect the measured intensity c(C), where £ 
is the etch depth, to depend on ry(z) and the initial NV 
concentration p(z) (before etching) according to, 



«t(C>0) = J dzr,(z-0p(z). (1) 



Thus, to obtain p(z) we must solve a deconvolution prob- 
lem. For our particular rj(z) a converging iterative solu- 
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tion can be obtained as follows: 



Po(z) 

Pn>o(z) 



^Pn{z) 



1 da 



'7(0) dC 
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(2) 
(3) 



dz 2 r)'(z 2 - z)p n - 1 (z 2 ) , (4) 



where rf{z) = drj/dz. The observed rate of convergence 
was such that, after 10 iterations, the remaining correc- 
tions had Max|/9„(z)|/Max|/9 (-z)| < 10~ 5 . In analyzing 
the data we must also contend with a large amount of 
noise. We believe this noise is due mainly to errors in 
replacing the sample before each measurement and in 
setting the focus, and we also have uncertainty in esti- 
mating the etch depth. The observed fluctuations are 
much larger than the Poisson detection limit and are ap- 
proximately 10% of the count rate. We used the fol- 
lowing procedure to estimate p(z): (1) perform a linear 
interpolation of log <j(z), and exponentiate to obtain a 
continuous version of p(z). (2) Calculate the correspond- 
ing p(z) using Eqs. [2]|4] (3) Smooth the result using a 
moving average with 100 nm total width to obtain p(z). 
(4) Calculate a ct(C) from p{z) using Eq. [I] and compare 
with the original data. The calculated p[z) are plotted in 
Fig. El and the corresponding tr(C) are the fits shown in 
Fig. [3] Including spontaneous emission modification had 
a noticeable effect on the calculated NV profiles, but the 
results are qualitatively similar to those obtained simply 
by differentiating a((). In order to estimate the absolute 
NV concentrations shown in the figure, we also made use 
of separate calibration measurements performed on sin- 
gle NV centers, where the count rate per unit excitation 
power was approximately 1.5 x 10 4 s~ 1 mW~ 1 , and the 
effective collection area of the confocal setup was deter- 
mined to be approximately 0.5 pm 2 . The estimated sheet 
densities of optically active NV centers for the three mea- 
surements are 1.3 x 10 13 cm -2 and 3.5 x 10 11 cm~ 2 for the 
high- and low-nitrogen sectors, respectively, of the HPHT 
sample, and 4.3 x 10 9 cm -2 for the CVD sample. 

The results in Fig. [4] confirm that, for the high-nitrogen 
sector of the HPHT sample, most of the optically ac- 
tive NV centers occur within the first 200 nm. The dis- 
tribution is significantly broader than the predicted va- 
cancy profile before annealing shown in Fig. [I] Such a 
spreading of the distribution might be expected because 
of vacancy diffusion during annealing. However, for the 
measurements on low-nitrogen material, the behavior is 
quite different. In these cases, the concentration of op- 
tically active NV centers is very low in the first 100 nm, 
and the maximum concentration occurs approximately 
250 nm below the surface. 
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FIG. 3: Integrated PL intensity vs. etch depth for the HPHT 
sample, high-nitrogen sector (red circles) and low-nitrogen 
sector (green squares) , and the CVD sample (blue diamonds) . 
For the CVD sample, the plotted intensity is the average in- 
tensity measured within the implanted squares minus the av- 
erage intensity outside of the implanted squares (see text). 
Curves: PL intensity computed from the NV concentration 
estimates shown in Fig. [4] 



IV. SPECTROSCOPY RESULTS 

Next, we present spectroscopy results from the HPHT 
sample that help explain the lack of PL originating near 
the surface in the low-nitrogen material. The insets of 
Fig. [5ja,b) show typical room-temperature PL spectra 
from the high- and low-nitrogen sectors. The most im- 
portant features are the 575 nm and 638 nm peaks which 
have been identified as the zero-phonon lines (ZPLs) of 
NV° and NV~ , respectively^. The sharp line at 572 nm 
that appears prominently in Fig. [5jb) is the diamond Ra- 
man line. To extract the NV~ and NV° PL intensities 
from each spectrum, we fitted a Voigt approximation^ 
to the 575 nm and 638 nm peaks after subtracting a lin- 
early sloped background fitted to regions on the two sides 
of the peak. The extracted peak areas are plotted as a 
function of etch depth in the main plots of Figs. [5ja,b). 
The NV°/NV~ intensity ratio for both sectors is plotted 
in Fig.^c). 

For the first few hundred nanometers of etching, quite 
different behavior is observed for the high-nitrogen and 
low-nitrogen sectors. In the high-nitrogen sector, the 
NV°/NV _ ratio is approximately 0.2 and changes lit- 
tle as the sample is etched, although it eventually de- 
creases. However, for the low-nitrogen sector, this ratio 
changes dramatically as the sample is etched, starting 
out at 2.2 after the first successful etch, and decreasing 
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FIG. 4: (a) Calculated Purcell factor (blue) and correspond- 
ing relative collection efficiency (red) for a [lll]-oriented NV 
center close to a (100) diamond surface, (b-d) NV concentra- 
tion vs. depth, estimated from the data in Fig. [3] 



to 0.3 after 400 nm of etching. Even more surprisingly, 
for the first 150 nm of etching the NV~ intensity actu- 
ally increases. The measured PL intensity includes the 
contributions from NV centers at all depths, and we do 
not expect significant absorption in the top layer. We 
believe the most likely explanation is that the top layer 
interacts electronically with the layers below. A simple 
model describing this interaction is presented in the next 
section. 

Finally, we briefly address the low-temperature spec- 
tral linewidth in the implanted samples, a critical issue 
for some potential applications. In the ideal case, the 
linewidth of NV centers created through ion implanta- 
tion and annealing would be equal to the linewidth ob- 
served in background NV centers distributed throughout 
the sample. For another Element 6 CVD sample with 
nitrogen content similar to that of the sample used in 
the etching experiment, the low-temperature linewidth 
(FWHM) of the NV" ZPL at 637 nm was measured to 
be 0.035 nm for the background NV centers. For the NV 
centers created by ion implanation (10 12 cm~ 2 , 200 keV 
Ga + ) followed by etching to remove the top 300 nm, the 
linewidth was 0.057 nm, and the tails of the ZPL peak 
were more prominent than for the background NV cen- 
ters. For the HPHT sample, the high-nitrogen sectors 
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FIG. 5: (a,b) Estimated NV° and NV" zero-phonon line 
(ZPL) intensities vs. etch depth for the (a) high-nitrogen 
and (b) low-nitrogen sectors of the HPHT sample. Insets: 
example spectra showing the diamond Raman line at 572 nm 
and NV° and NV - zero-phonon lines (ZPLs) at 575 nm and 
638 nm, respectively, (c) The NV°/NV~ intensity ratio vs. 
etch depth for the high- and low-nitrogen sectors of the Sum- 
itomo sample. 



show broad linewidths ranging from 0.5 — 1.9 nm. In one 
of the special low-nitrogen sectors, a strain-split NV - 
ZPL with 0.13 nm total linewidth was observed under 
normal excitation conditions. However, when collect- 
ing PL through a GaP waveguide using evanescent cou- 
pling 18 a much broader linewidth (1.0 nm) along with a 
large NV°/NV~ ratio was observed at the same location. 
This indicates that NV centers close to the surface had 
severely degraded properties. We have performed low- 
temperature spectroscopy on many samples implanted 
by a variety of methods and plan to report the results 
elsewhere^. 
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V. SIMPLE DEPLETION MODEL 

In this section we seek to explain the behavior in 
Fig. [j^b), especially the initial increase in the NV - sig- 
nal observed during the first 150 nm of etching. In our 
view the most likely explanation for this behavior is an 
electronic interaction, and in particular the presence of 
an electronic acceptor layer near the surface that removes 
electrons that would otherwise be available to form NV - 
in the layers underneath. 

It has been shown that the relative PL intensities of 
NV° and NV - depend on the concentrations of impuri- 
ties su ch as nitr ogen, and on the optical excitation con- 
ditioniPsHsMEHl Several types of impurities are present 
that can act as electronic donors and acceptors. Ni- 
trogen impurities are expected to be the main electron 
donor, with an energy level 1.7 eV below the conduction 
band minimum^. The NV center itself is known to ex- 
ist in the NV~ and NV° states, and perhaps could also 
exist in an NV + state. Because NV - has discrete ex- 
cited states 1.94 eV above the ground states, the ground 
states of NV - must be > 1.94 eV below the conduction 
band minimurrPSl. An abrupt change reported in the 
NV°/NV~ ratio at an excitation wavelength of 480 nrrPHl 
suggests that the actual position of the NV - ground 
state is 2.58 eV below the conduction band minimum. 
Finally, we can only explain the results in Fig. [5] if some 
species of electron acceptor is present in the top layer. 
As discussed below, there are several possible acceptor 
species. For the present discussion we label the accep- 
tor as A - , and choose an energy level 1.4 eV above the 
valence band maximum as expected either for substitu- 
tional Ga impurities^ or for graphitic defects^. The 
indirect bandgap of diamond is E g = E c — E v = 5.5 eV. 
Theoretically, for a junction between a strongly "p-type" 
layer and an "n-type" layer (nitrogen as the predomi- 
nant impurity), a "depletion layer" of ionized nitrogen 
impurities is expected to have a width of approximately 
x n = y // 2eAE/[N]e = 90 nm if the nitrogen concentra- 
tion is 1 ppm. Thus it is quite possible that a depletion 
layer could be responsible for a reduced NV - concentra- 
tion before etching in the low-nitrogen material. For a 
nitrogen concentration of 100 ppm, the theoretical deple- 
tion width is only 9 nm, and interactions between differ- 
ent layers would be difficult to observe. 

As a first step toward quantitative modeling of the 
spectroscopy results, we have implemented a simple 
model that solves the 1-D Poisson equation, 

= = f (n-[N+] + [NV-] + [A"]-p) , (5) 

under the condition that the electronic occupations of 
each type of impurity, as well as the conduction-band 
electron concentration n and valence-band hole concen- 
tration p, arc in thermal equilibrium according to Fcrmi- 
Dirac statistics. Eq. [5] is then solved numerically under 
the constraint that the electric field is zero at z = ±oo, or 
equivalently, that the total charge is zero. In the actual 



experiment, we do not necessarily expect the electronic 
populations to be in thermal equilibrium, since the sam- 
ple is subject to optical excitation at 532 nm. Neverthe- 
less, given our lack of knowledge about all of the optical 
excitation rates, this model is a good starting point. In 
order to simulate the experiment, we must also assume 
some distribution of nitrogen impurities, NV centers and 
acceptors. As an example, we have simulated the fol- 
lowing distribution, where z > is the depth below the 
original surface: 

[N] = [N°] + [N+] = N bkg - [NV] , (6) 

[NV] = [NV - ] + [NV ] = T^y N bkg , (7) 

r(z) = r o (l + 0(z-z o )(e- (z ~ Zo)/Az -l)) , (8) 
[A] = [A-] + [A ] = a 6( Zl - z) , (9) 

where Nbkg is the initial nitrogen concentration (assumed 
to be uniform), r(z) is a dimensionless rate for conver- 
sion of N to NV with a depth profile controlled by the 
parameters ro, zq and Az, 9{z) is the unit step function, 
and do is the concentration of acceptors in a top layer of 
thickness z\. We have chosen this impurity distribution 
because it is simple and yet contains the most essential 
features of the NV distribution and acceptor layer as dis- 
cussed above. 

The simulations shown in Fig. [6] used the parameters 
Nbkg = 1-76 x 10 16 cm- 3 (0.1 ppm), r = 15, z = Az = 
100 nm, z x = 150 nm, and a = 1.4 X 10 16 cm" 3 . The 
temperature was chosen to be 300 K, which produces es- 
sentially the same results as OK. Fig. |6ja) shows the 
calculated 2-D concentrations of NV - and NV° (the 3- 
D concentration integrated over z) as a function of etch 
depth. For these parameters, the simulation produces a 
similar behavior as observed in the spectroscopy data, 
where the NV - PL intensity increases at first, and then 
decreases. To help in explaining this behavior, energy 
band diagrams and impurity electron concentrations are 
shown before etching (b,c) and after etching 200 nm (d,e). 
Before etching, the 150 nm p-type layer takes electrons 
from nearby NV centers and nitrogen donors, so that NV 
centers cannot exist in their negatively charged state to 
a depth of 300 nm. Below this depth, NV - can exist, but 
the total NV concentration there is low. However, after 
removal of the p-type layer, nearly all of the remaining 
NV centers can exist in their negatively charged state. 

Thus, through a suitable choice of parameters, we can 
obtain good qualitative agreement with the spectroscopy 
results shown in Fig. [5jb) for the low-nitrogen sector, 
demonstrating that electronic depletion is a plausible ex- 
planation for the initial rise in the NV - signal. Such 
agreement is possible only for a limited range of param- 
eters. A substantial increase in NV - due to removal of 
the p-type layer is possible only if Nbkg is small, and 
Nbkg cannot be much larger than the value of 0.1 ppm 
used here. This in turn limits the total number of NV 
centers that can exist in the model. The maximum sheet 
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FIG. 6: Simulation of electron statistics using a thermal equi- 
librium model (see text for parameters): (a) Total NV° and 
NV~ populations (integrated over z) vs. etch depth, (b) Band 
diagram before etching, showing the conduction band mini- 
mum (c.b.), valence band maximum (v.b.), nitrogen donor 
level (N°), NV level (NV - ), and acceptor level (A - ), rela- 
tive to the Fermi level (at zero energy), (c) Populations be- 
fore etching of neutral nitrogen impurities (N°), negatively 
charged NV centers (NV - ), neutral NV centers (NV°), and 
acceptors (A) vs. depth z. (d) Band diagram after etching 
200 nm. (e) Populations after etching 200 nm. 



density of NV - obtained in the model, 3.4 x 10 11 cm -2 , 
is in good agreement with the experiment: as discussed 
above, we estimated 3.5 x 10 n cm -2 . In order for the 
NV - population to increase gradually as the p-type layer 
is etched, the acceptor concentration must also be small, 
comparable to the initial nitrogen concentration in this 
case. If it is much higher, the NV - population remains 
approximately constant until the p-type layer is almost 
completely removed, at which point the NV - popula- 
tion abruptly increases. The acceptor concentration used 
here, approximately 0.8 ppm in a 150 nm layer, corre- 
sponds to a sheet density of only 2.1 x 10 n cm -2 , or 
only 7% of the implanted Ga dose. 

One important difference between theory and experi- 
ment can be seen in the NV°/NV" ratio. In the model, 
this ratio falls rapidly to zero once the damaged layer 



has been removed, and the NV concentration in the re- 
maining layers is less than the nitrogen concentration. 
In the experiment, this ratio reaches a minimum of ap- 
proximately 0.25, and in general we observe this ratio to 
be dependent on the optical excitation power. This is 
a limitation of the thermal equilibrium model presented 
above, especially at low temperature. As the tempera- 
ture in the model is increased (a crude way to simulate 
intense optical excitation), the ratios will change, and the 
equilibrium value will eventually depend on the effective 
masses of electrons and holes in the conduction and va- 
lence bands. To model the system more accurately, we 
would need to use a rate-equation model for optical exci- 
tation combined with drift and diffusion of electrons and 
holes in the conduction and valance bands. Such a non- 
equilibrium model is feasible to solve numerically, but 
the main difficulty is that we have little knowledge of the 
ionization rates of nitrogen impurities, NV centers, and 
acceptors under optical excitation, as well as the electron 
and hole capture rates for these impurities, all of which 
are important. We have observed time-dependent behav- 
ior in the low-nitrogen sector that gives some hint of these 
rates: for excitation with 4.5 fiW at 532 nm we observed 
that if the excitation is suddenly switched on, the PL in- 
tensity is initially smaller and increases by a factor of two 
with ~ 70 ms rise-time. This rise time could be associ- 
ated, for example, with optically induced charge transfer 
from acceptors to NV centers in the depleted region by 
means of free carriers. 

Finally, let us briefly address the question of what ac- 
ceptor species could be responsible for creating such a de- 
pletion layer. Several possibilities exist. Damage close to 
the surface, created by polishing, for example, could have 
been present even before ion implantation. If such dam- 
age consists of graphitic defects, this could pin the Fermi 
level 1.4 eV above the valence band maximum 52 . How- 
ever, such defects would have to be distributed through- 
out the top ~ 150 nm of diamond, in a low concentra- 
tion, to explain the observed result. Ion implantation in- 
troduces vacancies, self-interstitials, and Ga impurities. 
The subsequent annealing process is complex, and the 
final end-products are uncertain. Monovacancies could 
remain, but with an acceptor level apparently positioned 
2.85 eV above the valence band maximurrPSl (nearly iden- 
tical to NV - ) , it is difficult to explain the observed results 
in terms of this acceptor. However, if vacancies combine 
to form extended defects^, a high concentration of accep- 
tor states throughout the diamond bandgap could exist. 
Finally, it is quite possible that the implanted Ga atoms 
are the dominant acceptor. Initially we expected that 
Ga impurities would not be important, since according 
to the simulation in Fig. [T] ~ 770 vacancies are produced 
for each implanted Ga + ion. However, the sheet density 
of acceptors needed to produce a depletion layer in our 
model is only 7% of the implanted Ga + dose. Accep- 
tor levels for substitutional Ga impurities and for Ga- 
vacancy complexes have recently been calculated, with 
positions 1.4 eV and 1.7 eV above the valence band max- 
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imum, respectively^. 

VI. CONCLUSIONS 

From the data and simulations presented above, we 
can draw several conclusions related to formation of NV 
centers through implantation of heavy ions followed by 
annealing. First, the vertical distance over which vacan- 
cies diffuse to form NV centers is only a few hundred 
nanometers, even for a fairly high annealing temperature 
of 925° C. Second, if the implantation dose is chosen 
to maximize the PL intensity in diamond with a mod- 
erately low nitrogen concentration (~ 1 ppm) , most of 
the NV~ PL that is produced originates from NV cen- 
ters > 200 nm below the surface, while NV centers closer 
to the surface are either in the NV° state or in other 



states that do not emit in our detection band. Such a 
situation is not useful for coupling NV centers to opti- 
cal structures. Finally, the spectroscopic results in the 
low-nitrogen diamond are highly suggestive of electronic 
depletion effects. Our model shows that the observed 
behavior can be explained by the presence of electronic 
acceptors near the surface, which could consist of impu- 
rities or defects introduced through the ion implantation 
and annealing, or damage originally present in the sam- 
ple. Controlling such effects will be important in finding 
improved methods to couple dense ensembles of NV cen- 
ters to optical structures on a diamond surface, and also 
for fabrication of single negatively-charged NV centers 
in ultra-pure diamond, where electronic depletion layers 
thicker than 1 /im are possible. 

The authors thank Dirk Englund and Steven Prawer 
for helpful discussions. 
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